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This paper describes a study on fractional flux periodicity of the ground state in planar 
systems made from a square lattice whose boundary is compacted into a torus. It is pointed 
out that, in the long length and large diameter limit of a torus, the ground-state energy 
and persistent currents show a fractional period of the fundamental unit of magnetic flux 
depending on the twist around the torus axis. We discuss the possible relationship between 
the twist and genus of a torus, and comment on fractional flux periodicity in a cylinder. 

§1. Introduction 

The fact that a single electron wave function has a fundamental unit of magnetic 
flux = 27r/e **' is known to be an essential consequence of quantum mechanics 
and is manifest in the Aharonov-Bohm effect. 1 ' The solid state properties of mate- 
rials are governed by many electrons, and each constituent has the <?o periodicity. 
However, the fundamental flux period of a material is not always <?o- F° r instance, 
superconductors exhibit a period of #o/2, which can be understood by charge dou- 
bling due to the Cooper pair formation. This must be a trivial fact if one suppose 
that a fundamental particle (or quasiparticle) has a — 2e charge due to attractive 
interaction, and that the fundamental flux period of a material is equal to that of 
the quasiparticle in the system (2ir/2e). 

In Fig. ^ we depict schematic diagrams of typical experimental settings for 
the Aharonov-Bohm effect. Consider an electron going over and under a very long 
impenetrable cylinder, denoted as a black circle in Fig. ^a). In the cylinder there 
is a magnetic field parallel to the cylinder axis, taken as normal to the plane of the 
figure. The probability of finding this electron in the interference region depends 
on the magnetic field ^ and the interference pattern having period ^o- When an 
electron is placed in a ring pierced by a magnetic flux ^ as shown in Fig. O^b), 
a current I defined by differentiating the ground state energy with respect to the 
magnetic flux, appears and is expected to have the flux periodicity ^o- 

Actual mesoscopic rings contain many electrons, and non vanishing currents, 
known as persistent currents, 2 ^ 4 ^ are also expected. Although the ground state 
in mesoscopic rings consists of many electrons, because each electron has the flux 
periodicity @q one might think the flux periodicity of the systems is the same as the 
flux periodicity of the constituent electron. In this paper, we present an example 
in which the fundamental flux period of the ground-state does not coincide with 



*' E-mail address: sasaken@imr.tohoku.ac.jp 
**' The electron charge is denoted by — e. We will use the units of h — c — 1. 
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Fig. 1. Schematic diagram of (a) an Aharonov-Bohm experiment and (b) Aharonov-Bohm ring. In 
(a), electrons are emitted from the source region, denoted by the empty circle, and produce the 
interference pattern. In (b), one- dimensional ring is pierced by a magnetic field and exhibits a 
persistent current. 



the flux periodicity of the constituent particle. We will show that a fractional flux 
periodicity (^q/2, <&o/3,- • •) can be realized in a certain limit of planar systems made 
from a square lattice whose geometry is a torus *) . To analyze the flux periodicity of 
the systems, we consider the ground state energy and persistent currents. These are 
considered suitable for examining the flux periodicity of the ground state since they 
are regarded as the Aharonov-Bohm effect in solid state systems. 4 ) In the previous 
paper, 5 ) we showed a typical example of a fractional periodicity in the Aharonov- 
Bohm effect. Here, we extend the theory in general cases and give detailed derivation 
of the results. 

This paper is organized as follows. In Sectional we obtain an explicit formula 
for the ground state energy and plot it for different lattice structures, showing flux 
periodicity numerically. In Section |*3 we examine the persistent currents and prove 
fractional flux periodicity in the long length and large diameter limit of a torus 
noting that the characteristic features of the ground state energy can be understood 
in terms of the twist-induced gauge field. In Section |1J we give a discussion and 
summary of our results. In Appendix we prove a formula used in the text and 
estimate a correction to a fractional periodicity. Throughout this paper we ignore 
electron spin and the effect of finite temperature, and consider a half-filling system. 

§2. Ground state energy 

The lattice structure of a torus is specified by two vectors: the chiral **' , given 
by Ch = NT X + MT y , and the translational, given by T w = PT X + QT y , where 
N, M, P and Q are integers and T a {a = x, y) is the unit vector of the square lattice 

*' It is noted that fractional (^o/3, 4>o/4, • • • except <2>o/2) flux periodicity presented in this 
paper holds in the long length and large diameter limit of a torus. 

**) We borrow this terminology from the carbon nanotube context. 6 ^ 
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in the direction of x and y where we define x and y in the two dimensional map 
(see Fig. Efb)). Here, Ch (T w ) characterizes the lattice structure around (along) the 
tube axis. We rewrite T a in terms of the chiral and translational vectors for a later 
convenience as 

'T x \ 1 / Q -M\ (C h s 



T y J N s \-P N ) \T W ' ' 1 ! 

where we define N s = NQ — MP, and |iV s | corresponds to the number of square 
lattices on the surface of the torus given by Ch x T w = N S (T X x T y ). In this paper, 
we take the x axis in the direction of Ch [M = 0), and denote P = 5N to describe 
the twist around the tube axis. Then, Ch • T w = 5NNa 2 where a(= \T X \ = \T y \) 
is the lattice constant. We depict an example of twisted torus in Fig. [21 in which 
all sites are connected by a line consisting nearest-neighbor bond. Fig. I^a) shows 



Fig. 2. Lattice structure of a twisted torus. The lattice configuration is given by N = 6, Q — 40 
and SN = 1. 

an illustration of the twisted torus. A torus can be unrolled to a parallelogram 
sheet as shown in Fig. EKb) . In the figure, the two lines extending upward from 'u' 
and downward from 'd' and having the same 'x' at the junction are not joined for 
a twisted torus, as shown in the inset of Fig. Eta). There are SN square lattices 
between the two lines. 

We consider the quantum mechanical behavior of conducting electrons on a 
square lattice using the nearest-neighbor tight-binding Hamiltonian: 

n = -t^E a l a e- ieAC ^a l + h.c, (2) 

i a=x,y 

where t is the hopping integral, A ex is a constant external gauge field, and aj and a\ 
are canonical annihilation-creation operators of the electrons at site i that satisfy the 
anti-commutation relation {aj,aj} = 6ij. We diagonalize TL in terms of the Bloch 

basis parametrized by wave vector k and obtain H. = ^2 k E(k — eA e *)a\ak, where 
afc = Y^i elk Tia i ^ s the annihilation operator of the Bloch orbital. The energy 

eigenvalue of the Hamiltonian E(k — eA cx ) is given by, 

E{k - eA cx ) = -2m ( ^ e^ k - eAex > Ta J . (3) 

\a=x,y ) 
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Fig. 3. Schematic diagram of (a) a twisted torus (upper) and the extension around a cross section 
(lower), and (b) its net diagram. It is convenient to consider a parallelogram sheet of square 
lattice as the net diagram of a twisted torus, where the chiral and translational vectors become 
the vectors at the sides of the parallelogram. The lattice configuration is given by N = 14, 
Q = 16 and 5N = -3. 



The Bloch wave vector satisfies the periodic boundary condition for Ch and T w 
through which the geometrical (or topological) information of a torus, such as the 
twist, is entered into the energy eigenvalue. We decompose wave vector k as fiiki + 
where k\ and ki are defined by 

Ch • k\ = 2tt, Ch • ki = 0, T w • k\ = and T w ■ k% = 2ir. (4) 

Here, k\ and ki can be rewritten in terms of the reciprocal lattice vectors K x and 
K y as 



fcl 



1 



(QK X - SNK : 



V)i 



A: 2 



1 



(NK y 



(5) 



where the reciprocal lattice vectors are defined by K a ■ Tb = 2n8 a b with a,b = x, y. 
Components [i\ and [ii are quantum numbers (integers) and define the Brillouin 
zone as 

~Q~ 
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N 
~2 



+ 1 < Mi < 



~N~ 
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+ 1 < M2 < 



(6) 



where [n] indicates the maximum integer smaller than n. 

We point out that the wave vector shifts according to the twist of a torus. For 
example, in the absence of twist (5N = 0), k\ and ki are proportional to K x and 
K y , respectively and perpendicular to each other; otherwise, they are not. The 
shift of the wave vector plays an important role when we consider the ground state 
configuration, which we discuss below. Here we rewrite Eq.(j3J) by means of Eqs.Q 
(with M = and P = 5N) and ® as 



E(k - eA ex ) 



-2t 



cos 



27T/i] 



+ cos 



27r( / u 2 



SN 
N 



Mi) - eA e 



Q 




(7) 



where we assume A ex ■ Ch = 0, and A eyi ■ T w corresponds to the Aharonov-Bohm 
flux penetrating the center of the torus. The ground state is defined as the state for 
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Fig. 4. Brillouin zone and Fermi lines in momentum space (a), and schematic diagram of the wave 
vector shifting due to the twist around the axis (b). There are two factors we should take into 
account to define the ground state. The first factor is the effect of the twist. Compared with 
the ground-state configuration in an untwisted torus, Eq.JSJ indicates that the location of the 
energy eigenstates shifts depending on twist SN and on the electron motion around the axis 
(fii). This shift is illustrated in (b), where we denote the Bloch states as circles and depict 
several states of positive negative hi and /ii = 0. In the absence of twist, all states are 
located at mesh line crossings. When twist is present, the states shift depending on their \x\ 
and SN/N. The second factor is the external gauge field. The external gauge field can shift the 
wave vector, however, the shift is independent of electron motion contrary to the case where 
twist is present. 



which all electronic states below (above) the Fermi surface are occupied (empty). We 
consider a half-filling system and set the Fermi energy E-p = 0. We define the Fermi 
line in the two dimensional Brillouin zone, within which all states are occupied, as 
the solution of E(k) = 0, forming a square as shown in Fig. Ufa). The side of the 
square (or Fermi line) is denoted by /c ± and k±. Each energy eigenstate is labeled 
by fj,\ and fi2 and thus we can specify the occupied states by /Ui and fi2 explicitly. 
In fact the region inside the Fermi line is defined as 

. 27r/ii ^ 2vr( j » 2 -^ / i 1 )-eA ex -T w 

+ q (8) 

. 2ttmi 2^(,» 2 - ffi/xQ - eA™ ■ T w 
"-< — ^ <vr, (9) 

where we note that the Fermi lines satisfy ■ (T x +T y ) = ±7r and k± ■ (T x — T y ) = ±7r, 
and Eqs.(jHJ) and ® can be rewritten as 

kT ■ (T x + T y ) <{k- eA cx ) ■ (T x + T y ) < k + ■ (T x + T y ), (10) 
k. ■ (T x - T y ) < (k - eA ex ) ■ (T x - T y ) < k + ■ (T x - T y ). (11) 
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For a fixed value of we obtain the following region of M2 ; in which all states are 
occupied by electrons: 



^ { 2 N 



for m > 0, and 



" 1 2 N 



+ —/xi + iV* 



+ —Ml + N 



+ 1 < M2 < 



+ 1 < M2 < 



1 , Mi 

2 iV 



+ —mi + iV* 



(12) 



(13) 



for mi < 0, respectively. Here, we set A ex ■ T w = and define the number of 
Aharonov-Bohm flux N& = ^/^o- The ground-state energy is therefore 



[f] 

E 



1 [Mil 

2 2V 



" AT 



E 



- e,4 ex ). 



(14) 



This is the explicit mathematical expression for the ground state energy of a twisted 
torus. 
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Fig. 5. Ground state energy as a function of flux. We fix TV = 80 and Q = 160, then plot for 
SN — 1,20,40 to examine the dependence of the ground-state energy on twist. The observed 
flux periodicity seems to be #o/2 for SN — 1 (solid line) and SN = 40 (dashed line) and <2>o/4 
for SN = 20 (dotted line). For the untwisted torus (N = 40, Q = 160 and SN = 0), refer to 
Fig- Egl 



Here, we plot Eo(^)/At for different N, SN and Q as a function of N&, and 
observe the characteristics of the results. Fig. El shows the dependence of the ground 
state energy on SN, where we fix N = 80 and Q = 160 and vary the twist as 
SN = 1,20,40. The flux periodicity seems to be <P /2 for SN = 1 (solid line) 
and SN = 40 (dashed line), and <£o/4 for SN = 20 (dotted line). It appears that 
flux periodicity of the dashed and dotted lines corresponds to the lattice structure of 
5N/N. We note that ^o/2 periodicity of solid line (5N = 1: odd number) is accurate 
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in our numerical data (10-figure numbers) and that of dashed line (5N = 40: even 
number) is approximate one. In addition to observing flux periodicity, we see that 
the solid line is smooth (no cusp) at N& = for example, when compared with the 
other two lines. 

Next we examine the dependence of the ground state energy on the diameter of 

6(c) and |6(d)| show the ground state energy with 
20 while varying N = 40,80,81, 160. Observing the flux 



the tube N. Figs. |6(a)[ |6(b) 

fixed Q 



160 and 5N 

periodicity of the samples, they appear to be periodic in Figs. 6(a) and |6(b)| with a 
period of <&o/2 in Fig. 6(a) and ^o/4 in Fig. 6(b) They are not periodic, however, 
due to a small modulation in the ground state energy. This is manifest in Figs. |6(c 



and |6(d)| This modulation can be thought of as a finite size correction which will 
be discussed in Section |3J 
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(a) Ground state energy for N = 40, Q = 
160 and SN = 20. We plot (E ($) - 
E (0))/4t where E {0)/4t = -1296.892. 



(b) Ground state energy for N = 80, Q 
160 and SN = 20. We plot (E (&) 
E (0))/U where E (0)/U = -2593.812. 
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(c) Ground state energy for iV = 81, Q = 
160 and SN = 20. We plot (E ($) - 
E {0))/4t where E (0)/4t = -2626.2443. 



(d) Ground state energy for N = 160, 
Q = 160 and SN = 20. We plot {E {$) - 
E (0))/U where E (0)/4t = -5187.640. 



Fig. 6. Dependence of the ground state energy on the diameter of the tube N. 



Finally, we observe the dependence of the ground state energy on Q shown in 
Figs [TfaJ} [7(b)| [7(c)] and [7(d)| Here, we note that all previous examples except 
Fig. 6(c)| have had the common feature of Q as a multiple of N. Fig. 7(a) shows 



*' This figure is a refined version of the dotted line in Fig. |K| 
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the standard flux periodicity <Pq. In Fig. 7(c)[ however, we have an approximate 



periodicity of &o/2 even in the absence of twist. We note that &0/2 periodicity in 
Fig. [7(d)1 (Q = 181: odd number) is accurate in our numerical data. 




(a) Ground state energy for TV = 40, 
Q = 160 and SN = 0. We plot (E (&) - 
E Q (0))/U where E Q (0)/4t = -1296.83. 



^ 0.005 




(c) Ground state energy for N = 40, 
Q = 180 and SN = 0. We plot (E ($) - 
E (0))/U where £ (0)/4t = -1459.008. 



^ 0.002 




(b) Ground state energy for N — 40, 
Q = 170 and SN = 0. We plot {E ($) - 
E {0))/U where E {0)/4t = -1377.965. 




(d) Ground state energy for N = 40, 
Q = 181 and SN = 0. We plot (E ($) - 
E (0))/4t where E Q {0)/4t = -1467.135. 



Fig. 7. Dependence of the ground state energy on the diameter of the torus Q. 

It is important to mention that flux periodicity and the amplitude of the oscil- 
lation of Eq{$) depends much on N, Q and SN. We will discuss the dependence of 
the flux periodicity and amplitude on them in the next section. 

§3. Persistent currents and interpretation of the numerical results 

In the previous section, we presented a formula for the ground state energy. By 
plotting the energy for different values of N, 8N and Q, we found that some tori 
exhibit an approximate fractional flux periodicity. In this section, we will present an 
analytical results that demonstrate and prove fractional flux periodicity in the long 
length (Q S> 1) and large diameter (N 3> 1) limit for a twisted torus in terms of 
the persistent currents. By so doing, we clarify the relationship between the lattice 
structure and flux periodicity. 
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Persistent current is defined by differentiating the ground state energy with 
respect to the magnetic flux <£: 4 ) 

Because Eq(<P) consists of many energy dispersion relations as specified by \x\ (here- 
after referred to as \i\-Xh. energy band), we can write 

[f] 

E (<P) = E (Pl>*) ( 16 ) 

where E(/j,i,$) is the ground state energy of ni-th. energy band and is the sum of 
the energy at /X2-th k points which satisfy Eas.(|12|) and 

S( M i,*)= Eimh + ^h-eA™). (17) 

w=[-o(i- J ^ i )+^w+w*]+i 

Hence, the persistent current can also be decomposed into the current of each energy 
band as 

[#] 

Jpc(#)= /(/.!,#), where 1(^,3) = - gg^l*) , ( i 8 ) 

W =[-f]+l 

To calculate the persistent currents in long systems with |T W | 3> a (or Q 3> 1), it 
is not necessary to sum over the energy eigenvalue of the valence electrons but just 
calculate the Fermi velocity of the energy band. Because, when we vary the magnetic 
field, the change of the ground state energy is almost determined by the change of the 
energy eigenvalue of electrons near the Fermi level and that of the energy eigenvalues 
deep in energy bands cancel to one another. Using Fermi velocity, the amplitude of 
I(/j,i,<I>) (the persistent current for the ni-th. energy band) is well approximated by 4 ) 

iM-Zgp, <i») 

where \T W \ is the system length and Vp(fj,i) denotes the Fermi velocity of the (ii-th 
energy band. The proof of this approximation is given in Appendix^] By fixing \i\ 
and expanding Eq.Q on \ii around the Fermi level, we obtain the energy dispersion 
relation of the /ii-th energy band as 

W w = M»i)P2 + o( v l). (20) 



where P2(= ^2^2) denotes the momentum along the axis. Here, we neglect an exter- 
nal gauge field. Notice that Ea.(|20jl represents the case for [12 > 0, and for ^2 < 0, 
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the energy dispersion relation is given by 
the Fermi velocity and is obtained as 



-TL^ X up to <D(jp2). The coefficient of P2 is 



2ta 



sin 



2"7T/ii 

N 



(21) 



First we consider persistent currents in an untwisted torus (8N = 0). In this 
case, all energy bands produce the same function for the persistent current, which is 
a saw-tooth curve as a function of with the same zero points (the flux for which 
the amplitude of the current vanishes) but with different amplitudes. In fact, the 
amplitude of the total current is given by a summation of all amplitudes: 



'tot 



[-f]+! 



' 2eta 




7T 




cot 




IT 1 
, \ J -w | 
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eta 




IT 




cot 




IT 1 

V \ J -w I 




2N 



even, 



for N = odd. 



(22) 



Here, we implicitly assume Q is a multiple of N (see Appendix 0). The other cases 
will be mentioned in the end of this Section. Persistent current I pc in the torus is 
given by /tot multiplied by —(ft/ir ((ft = 2ir(<P/<I>o)) *) The linear relation between <P 
and the persistent currents has a periodicity of <Pq, and we then have a saw-tooth 
current: 



/to^B-ir 1 ^^. 



(23) 



n=l 



It is noted that fundamental periodicity obtained here is not an approximate 
periodicity but is exact, and also that persistent currents in a torus contain always 
<Po periodicity. However, the <Pq periodicity is not always a fundamental periodicity 
as is shown below. 

We comment on the functional form of the persistent currents. The saw-tooth 
persistent current is a characteristic of the non-interacting theories at zero temper- 
ature. Each saw-tooth curve loses its sharpness due to disorder, or at finite tem- 
perature. 7 - 1 We note here, and in all subsequent discussions, that we are implicitly 
assuming the amplitude of the persistent current vanishes at (ft = 0. Generally, the 
zero amplitude points depend on the total number of electrons in the system and 
may be (ft = 2irj or (ft = 2wj + it where j is an integer (Appendix ^J). As for the 
flux periodicity, the zero amplitude position does not affect the results. We note 
also that the persistent current corresponding to the ground state energy shown in 
Fig. 7(a) can be approximated by Ea.(|23|) with the replacement (ft —>(/> — vr. 

Next, we consider a twisted torus (5N ^ 0). To understand the effect of twist 
on the persistent currents, we introduce twist-induced gauge field ^4 twlst . In terms of 
^twist ^ seconc [ term of the right hand side of Eq. (J7J) can be rewritten as 



—2t cos 



2vr^ 2 - (eA cx + ^A twist ) ■ T u 

Q 



(24) 



*' We may multiply /tot by (fr/ir as well. The sign (±</>) depends on the choice of positive 
direction along the torus axis and does not concern flux periodicity. 
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can 



where we define ^ twlst ■ C h = and A twlst • T w = 2tt5N/N. Like A ex , A 
shift the wave vector. However, the shift depends on the twist (5N/N) and the wave 
vector around the tube axis (niki). Because of this, coupling between ^4 twlst and the 
conducting electron preserves the time-reversal symmetry of the whole system. This 
is contrasted with the time-reversal symmetry broken by A ex . The twist behaves as 
an extra gauge field and shifts the wave vector along the axis direction, so that the 
zero points of the persistent current (or a saw-tooth curve) also shift. 8 - 1 As a result, 
in order to calculate the total current, we must sum the saw-tooth curves that have 
different zero points and different amplitudes, all of which depend on \x\. Then we 
have 



2eta Ij N { 



(25) 



where iffi {(f) is expressed as 



WW = 2 - y f ( _i r+1 !i*t^ 

7r z — ' ^— ' n 

= I Y (-l)^i^M (l + cos (n,5N)) ' 

IT < ^ 71 



27T/ii 

2sin(f) 



n=l 



COS 



2irn&N \ 
TV ) 



COS 



» TV , 



, (26) 



where we assume N is an even number. When N is an odd number, we should add 
the following extra term to the right-hand side of Eq. (|26j) : 



^ oo 

a 

n=l 



1 



,n+l 



sin(n<j 



■ cos(nir5N) 



2sin(f) 



cos 



TV 



+ cos IfJ 



(27) 



Depending on an even or odd number of 5N, term 1 + cos(irn5N) in Eq. ()26|) vanishes 
if n5N is an odd number or is equal to 2 for other cases. Then Eo.(|26|) reduces to 



WW 



1 sin(2re^) co 
it ^ n n 

n=l 



2 E ( - 



1 



<n+l 



sin(n<j 



for SN = odd, 



-Ct for 5N = even, 



(28) 



7T 



n=l 



n 



where C° and are given 

2sin(f) 



Co 



COS 



' iirnSN - 
v N , 



COS 



2sin(f) 



cos 



'27m5TV> 
«. TV ) 



COS 



(29) 



The difference in periodicity for (ft is clear because sin(2ra<f;) has a period one half 
the fundamental unit of magnetic flux, or <Pq/2. We note that the one-half flux 
periodicity for an odd number of 5N becomes exact if N is an even number, which 
means that the ground-state energy Ea, H14j) has one-half flux periodicity (Q is not 
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necessary to be a large number). The solid line (SN = 1) in Fig. El is an example of 
5N = odd showing an exact one-half flux periodicity. This conclusion is valid for any 
even number of N. For an odd number of N, exact <Po/2 periodicity is absent when 
TV is a finite number because of Ea.(|27|). This accurate one-half periodicity is due 
to the fact that the persistent current is an odd function (see Appendix in the 
absence of twist, which provides the term 1 + cos(irn5N) in the persistent current 
even for a finite Q value. 

A fundamental flux periodicity for a current demonstrates a more interesting 
behavior for a specific combination of N and 8N. For a while, let us fix N an even 
number. We consider a particular structure for 5N/N = 1/3 where 5N is an even 
number. In this case, the fundamental period of is 3, i.e., C^ +3 = C^, which shows 
there may exist <&o/3 fundamental period in the currents. In fact, when N S> 1, we 
have Cf = Cf *J -8vr/3iV and Cf « 2iY/vr. We then obtain 



i s N N m = -YY sin((3rt + m)0) C* m 

n=u m=l 

=^E(- i ) n+isin( r 0)c7 3 e+0(i/iv) 

n=l 

2N 

« for iV > 1, (30) 

7T Z 



where in the last line of Eq . (|30j) we assume —tt/3 < <p < tt/3. The fundamental flux 
period of Eq. (j5U|) becomes <?o/3- It should be noted that this argument can be applied 
to other twisted structures. For instance, when 5N/N = 1/Z (1/Z = 1/4, 1/5, • • • ), 
we obtain another fractional period of We note that the total persistent 

currents are still saw-tooth curves as expected from the non-interacting theories. 
The ground state energy corresponding to Z = 2 and Z = 4 is plotted in Fig. El 
as dashed and dotted lines, where there remains a finite energy modulation which 
prevents an exact fractional periodicity. It is noted that, in the large diameter limit 
(N ^> 1), the same behavior holds for an odd number of N. 

Finally, let us remark that the amplitude of the current exhibits a nontrivial 
dependence on <P when N — > 00 with a fixed value of 5N. If we divide iffi ((f)) by N 
and take the limit of N — > 00, we then have *' 



lim 



7T 



2 sin(2n0) 



1 



71=1 

OO 



n 



1 - 4n 2 5N 2 



for 5N = odd, 



7T 2 ^ { 1 n 1 - n 2 5N 2 

n=l 



(31) 



even. 



*' In this case, Eq. l27t can be neglected because the correction to Eq. 13111 is given by 0(1/N 2 ). 
Therefore, the following argument can be applied regardless of an even and odd number of N. 
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When 5N S> 1, we sum n in the above equations and obtain 
r(vr-2</»)(7r 2 -(20-7r) 2 ) 



rSN 

lim 



(0 



7T 



2 24vr 2 2 5iV 2 



for < < 7T and &/V = odd, 



3tt 2 5N 2 



for — 7r < 6 < tt and <5iV 



even, 



where we have used the following mathematical formula: 
sin(n^) (vr - 0) (vr 2 -(</>- 7r) 2 ) 



E 

71=1 



??/' 
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n=l 



1 



,n+l 



sinlnc 



IV 



vr 



12 



for < <r5 < 2?r, 



for — TT < <j> < IT. 



(32) 

(33) 
(34) 



In Fig. El we plot the persistent currents of Eq.(j32J) for an even and odd number 
of 5N as a function of 4>. It should be noted that the persistent currents are not 
standard saw-tooth curves even though we are considering non-interacting electrons, 
and the functional shape of Ea, ()32|) for 5N is an odd number similar to sin(2</>) but 
not identical (see Fig. EJ). The solid line (5N = 1 = odd) in Fig. El may be regarded 
as an (approximate) example of Eq. (|3*Tj) and shows a smooth curve as a function of 
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Fig. 8. Persistent currents as N 

number of 5N. We multiply 2 (16) for an even (odd) number case 
(dotted line) for comparison. 



We also plot 0.8sin(2^>) 



summarizing our results *): (1) If N is an even number, there exists exact one- 
half flux periodicity when 5N is an odd number; (2) In the long length {Q ^> 1) 
and large diameter (N S> 1) limit, a fundamental flux period of the ground-state 
becomes fractional as $q/Z, depending on the ratio of twist SN to N (1/Z = 5N/N), 

*' These results are valid when Q is a multiple of N. 
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although we do not assume any interaction that can form a quasiparticle of charge 
Ze\ and (3) As N — > oo with fixed value SN, the currents are not standard saw-tooth 
as is normally expected from non-interacting theories at zero temperature. 

The primary factor in these results is the Fermi surface structure of the square 
lattice. It appears that many energy bands cross the Fermi level when we roll a 
planar sheet into a cylinder, and electrons near the Fermi level in each energy band 
contribute to the persistent currents which interfere with one another due to the 
twist. This point can be made clear by comparing a square lattice with a honeycomb 
lattice (which possesses only two distinct Fermi points). In this case, we can observe 
at most a <&o/2 periodicity. 8 ) 

Finally we briefly discuss the dependence of persistent currents on the value of Q. 
So far we implicitly assumed that Q is a multiple of N. For a general Q and a general 
structure (different chiral and translational vectors), the persistent currents are not 
as simple 8 ) as the results obtained up to now. For example, when the remainder 
of Q/N is an odd number, the one-half periodicity may emerge even in the absence 
of twist (see Fig. |7(d)| ). Furthermore, another fractional period may appear for 
a specific value of Q in a certain limit. This complication can be understood by 
Eas.(|12|l and (|13|). For /ii-th energy band, we have the occupied states 



+ 1 < M2 < 



(35) 



Suppose Q is an even number and we set Q = nN + 5Q (n is an integer). Then the 
above inequality reduces to 



Q 



HidN + \m\5Q 



N 



+ 1 <M2 



Q 



< I -j- -n\m\ ) + 



HiSN — \ni\6Q 
N 



+ N$ 



(36) 



which shows an asymmetry of the left and right Fermi point in the /xi-th energy 
band because of the (SQ/N)\/ii\ term. In fact, this term seems not to be regarded 
as the gauge field presented in this paper. However, if 5N = we can obtain several 
consequences from the point of view of twist-induced gauge field. We consider both 
Hi-th energy band in the first (second) quadrant and — /ii-th energy band in the 
third (fourth) quadrant of the Brillouin zone (see Fig.EJ), an d examine the persistent 
current contributed from those regions. In this case, for a fixed value of jX\, we obtain 
the following region of in which all states are occupied by electrons: 



— - n\m\ + 



+ 1 ■ (f 



n\n\\ + 



(37) 



This shows that the remainder 5Q can be regarded as a twist-induced gauge field and 
have a similar effect as the twist, which is manifest when we compare Eq.JHll) with 
Ea. ()35|) . Therefore, we can apply the same analysis to a non- vanishing remainder of 
Q/N even in the absence of twist. 
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§4. Discussion and summary 



Here, we discuss some possible extensions of our results. Because we have ob- 
served a fractional periodicity (in a certain limit) in the ground state, it may be 
interesting to ask "Is it possible that the persistent currents exhibit multiple periods 
such as 2<Pq or 3<Po as a fundamental period?" To answer this question, we consider 
higher genus materials (g: number of holes), the ground state of which exhibit g<Po 
periodicity depending on the genus (g = 2, 3, • • • ). 9 ^ Moreover, there are other exam- 
ples where the charge of a quasiparticle itself becomes fractional. The quasiparticle 
in the quantum Hall effect has fractional charge e/3 10 ) and there is a model con- 
taining a fractional e/2 charged soliton in 1+1 dimensions. 11 ) Those systems might 
exhibit 3<Po or 2<?o periodicity in the ground state, respectively. Considering the 
possibility that the system exhibits flux periodicity such as (2/3)^o> (2/5)^o ; when 
we consider a torus with g = 1 such a periodicity seems impossible. However, for 
higher genus materials {g > 2), there is still a chance for a fundamental flux period 
of (g/Z)<Po, where Z is an integer. Such a periodicity may exist when the twist and 
genus can be simultaneously defined for a higher genus material, although we have 
no idea if one can define twist in higher genus materials. 

We would like to mention the relationship between our results and the previously 
published literature on a nontrivial flux periodicity of persistent current. Cheung et 
al. 12 ) examined persistent currents in a finite height cylinder whose surface is made 
from two-dimensional square lattice. They found that when the cylinder is threaded 
by a magnetic flux, the persistent current is approximated by 



where Iq is a constant amplitude and M is the number of lattice site in the height 
direction and N in the circumference (see Fig. E[b)). They numerically checked for 
several samples that flux periodicity becomes (approximately) a fraction depending 
on the aspect ratio of height and circumference. It is valuable to note that Ea.([38|) 
becomes almost identical to Eg . 1)26(1 if we replace N — > 5N and 2(M + 1) = N. Thus, 
fractional flux periodicity can be verified when iV/2(M + l) = 1/Z in the large height 
M 3> 1 limit. The same configuration with a magnetic field applied perpendicular 
to the cylindrical surface was analyzed by Choi and Yi. 13 ) They reported that 
a fractional flux periodicity appears in the persistent currents depending on the 
perpendicular magnetic field and the number of lattice sites along the circumference. 
We also note that periodicity is observed in the conductivity 14 ) (not persistent 
currents) for cylindrical geometries. 

In this paper, we presented an example in which the fundamental flux period of 
the ground-state does not coincide with the flux periodicity of the constituent parti- 
cle. We showed that a fractional flux periodicity (^o/3, ^o/4, • • • ) can be realized in 
the long length {Q 3> 1) and large diameter [N 3> 1) limit of twisted tori made of 
square lattice for half- filling cases, and the correction to the amplitude of persistent 
currents scales as 0(N/Q 5 ). As for one-half flux periodicity (#o/2), we found that 




(38) 
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(a) $ (b) 

A M 




Fig. 9. Schematic diagram of (a) a twisted torus studied in this paper and (b) a cylindrical geometry 
studied by Cheung et al. 12 ' Both surfaces are made from a square lattice. 

if N is an even number, an odd number of 5N or 5Q gives an exact ^o/2 periodicity 
in the persistent currents. Table [H summarizes the relationship between the lattice 
structure and flux periodicity of persistent currents, derived in this study. First, we 
classify the tori into two types depending on whether or not 5Q is zero, where 5Q is 
defined as the remainder of Q/N. As we have discussed in the text, for SQ 7^ and 
5N = 0, we can obtain the similar results for 5Q = by setting 5N = 5Q. Further- 
more, for the case of iV —* oo with a fixed value of 8N{^$> 1), the persistent currents 
are not standard saw-tooth curves as expected from the non-interacting theories at 
zero temperature. 



Lattice Structure and Flux Periodicity 



SQ 


SN 


Period (in units of $o) 


SQ = 


SN = odd 


1/2 for TV = even (exact) 




SN 


1/Z for SN/N = l/Z (approximate) 



SQ <-> SN 



Table I. Nontrivial flux periodicity of tori. For an even number of TV, there exists one-half flux 
periodicity when <5TV (SQ = 0) or SQ (SN = 0) is an odd number. In the long length (Q 1) 
and large diameter (TV S> 1) limit, we obtain a fractional flux periodicity depending on SN/N 
or SQ/N. 
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Appendix A 

Persistent current and Fermi velocity 

We prove that the amplitude of the persistent current is approximated by Eq. (|19|) 
for a long system (Q 3> 1). For simplicity, let us first consider an untwisted torus 
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and fix Q as an even numbered multiple of N *'. Then 



^ l 2 N 



(A-l) 



holds for < Af<p < 1 because the right hand side of Ea. HA-l|) is an integer. It follows 
from Eqs.©,©,© and (lA~ll that for < N$ < 1 



Ml 



/ — 7T sm 



i+i 



2t 2vr /2vr|^i| 



sm 



27r(/i 2 - Afo) 
. - 2iV # ) 



sm 



^(l-2iV,) + (i, 



(A-2) 



(A3) 



where we have used Ea. (|21|) in the last line. This shows that the amplitude of the 
persistent current of /Ui-th energy band is approximated by Ea. (|19|) and a correction 
on the order of I(fii)0(l/Q 2 ) can be ignored for long systems of Q S> 1. Similarly, 
for an odd number of Q up to 0(1/Q) we have 



^-(-2N ) for -- < N$ < -. 



(A-4) 



Here, it is important to mention the functional form of the persistent current. The 
persistent current for the ni-th energy band, which is valid for any N$, is obtained 
by extending sin (|>(1 — 2N$)\ in J(/xi,#) fEa. (|A-2[l ) as an odd function of N&. 
Therefore, it can be expanded by the Fourier series: 



C n sin(27rniV$), 



(A5) 



71=1 



with an appropriate Fourier coefficient C n . This fact will be useful to understand an 
exact one-half flux periodicity, which is valid for an odd number of 5N (SQ) with an 
even number of N (see Eq. ()26(l and discussion below Ea.([29|)). 

Next, we consider a twisted torus where Q is an even numbered multiple of N. 
Then 



-Q 



iMil 
N 



-Q 



[Mil 



holds for < jffjLi + N$ < 1 because the right hand side of Eq.(|A 



(A-6) 

is an integer. It 



Note that throughout this Appendix, Q is fixed as a multiple of N. 



18 



follows from Ens.lf7jl.(fT7 ]l .ip |l and O that for (6N/N)$ < N$ < 1 - (6N/N)$ 

L J (^) = L L ^g sm q 



sin 



This shows that the amplitude of the persistent current of ^i-th energy band is ap- 
proximated by Ea. (fTT?|) and a correction on the order of I(/j,i)0(1/Q 2 ) can be ignored 
for long systems of Q 3> 1. As shown in Eas. (|23|) . (|3Uj) and (|31|) the summation of 
//i results in, roughly speaking, the change from siii(2tthi/N) to N. Hence, the 
correction to the persistent currents (I^f (4>)) is given by 0(N/Q 3 ). 
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